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ABSTRACT 

We present an observational study about the effects of the interactions in the kinematics, stel- 
lar population and abundances of the components of the galaxy pair AM 2306-721. Rotation 
curves for the main and companion galaxies were obtained, showing a deprojected velocity 
amplitude of 175 km s _I and 185 km s _1 , respectively. The interaction between the main and 
companion galaxies was modeled using numerical N-body/hydrodynamical simulations, with 
the result indicating that the current stage of the merger would be about 250 Myr after peri- 
galacticum. The spatial variation in the distribution of the stellar population components in 
both galaxies was analysed by fitting combinations of stellar population models of different 
age groups. The central region of main galaxy is dominated by an old (5-10Gyr) population, 
while significant contributions from a young (200 Myr) and intermediate (1 Gyr) components 
are found in the disk, being enhanced in the direction of the tidal features. The stellar popula- 
tion of the companion galaxy is overall much younger, being dominated by components with 
1 Gyr or less, quite widely spread over the whole disk. Spatial profiles of the oxygen abun- 
dance were obtained from the a grid of photoionization models using the /?23 line ratio. The 
disk of the main galaxy shows a clear radial gradient, while the companion galaxy presents 
an oxygen abundance relatively homogeneous across the disk. The absence of an abundance 
gradient in the secondary galaxy is interpreted in terms of mixing by gas flows from the outer 
parts to the center of the galaxy due to the gravitational interaction with the more massive 
primary. 

Key words: galaxies: general - galaxies: stellar content - galaxies: abundances - galaxies: 
interactions - galaxies: kinematics and dynamics - galaxies: starburst 



1 INTRODUCTION 

It is widely accepted by some time that merging and interaction 
events play an important role in the formation and evolution of 
galaxies. Mergers change the mass function of galaxies, creating a 
progression from small galaxies to larger ones; the merging process 
can also change the morphology of the constituents, transforming 
gas-rich spirals in quiescent ellipticals. Interactions can also trig- 
ger a wide set of physical and morphological phenomena, such as 
tidal tails, bridges and shells, kinematically decoupled cores and 
star formation enhancements (see review of lStruckll 1999T) . 

Interacting galaxies also show enhanced star formation 
when compared with isolated galaxies. S uch e nhancement was 
initially proposed by lLarson & Tinslevl (l978) to explain the 
wider range of optical colors found in galaxies in pairs. 
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Since then, numerous studies have confirmed these results, es- 
pecially in the central regions, through meas u rements of op- 
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sion ( Hummellll98ll) . Recent studies have also shown that this 
enhanceme nt is a function o f the projected galaxy pair sepa- 
ration (e.g. lBartonetaf]|2000l : lLambas et all 120031 ; iNikolic et ail 
|2004) , being stronger in low-mass than i n hig h-mass galaxies 
(e.g. IWoods & Gelled 120071 : lEllison et all |2008|) . In particular, 
Elliso n et alj J2008I) found that the star formation rate as measured 
by the Ho- equivalent width for galaxies in pairs selected from the 
Sloan Digital Sky Survey, is some 70% higher when compared to a 
control sample of galaxies with equal stellar mass distribution. 

There is a connection between the interaction strength and 
the morphological distortion in binary galaxies. According to 



Mih os & Hernquistl dl996h models, the response of the gas to a 
close passage depends dramatically on the mass distribution of 
the galaxy, with the irregularities in the gas velocity field trac- 
ing the disturbances in the gravitational potential of the galaxy 
as observed, for e xample, in some galaxies in the Virgo cluster 
jRubinetalJl 19991) . Combined N-body/hydrodynamic simulations 
show that galaxy-galaxy mergers disturb the gas velocity field sig- 
nificantly, an d hence lead to asymme tries and distortions in the ro- 
tation curves dKronberger et alj2006T) . However, according to these 
authors, no severe distortions are observable about 1 Gyr after the 
first encounter. 

The gas motions created by the interac tion can also signifi- 
cantly alter the ch e mical state of the galaxies dKoeppen & Arimotol 
Il990l : iDalcantonl |2007|) . and modify the usually smooth ra- 
dial metallicity gradient often foun d in isolated disk g alaxies 
faenrv & Worthevl Il999h . Recently, iKewlev et all ( 120061) found 
that O/H abundance in the central region of nearby galaxy pairs 
is systematically lower that that of isolated objects. These authors 
suggest that the lower metallicity is a consequence of gas infall 
caused by the interaction, but very few observational studies have 
been published analysing in detail the influence of different levels 
of interactions in the metallicity distribution and enrichment prop- 
erties of galaxies. 

Ferreiro & Pastoriza 2004 (hereafter FP04), in a study of the 
integrated photometry and star formation activity in a sample of in- 
teract ing systems from the Arp-Madore catalogue dArp & Madord 
fl987h . found that the galaxies involved have bluer colours than 
those of isolated galaxies of the same morphological type, indicat- 
ing an enhancement of star for ming activity. This enhancem ent was 
also previously suggested by Donzelli & Pastorizal d 19971) to ex- 
plain the slightly larger values of Ho- + N n equivalent widths found 
in these systems, when compared with normal isolated spiral galax- 
ies. From their sample, we have selected several systems to start a 
more comprehensive study of the effects of the interactions in the 
kinematics, stellar population and abundances of the galaxies in the 
so-called "minor merger" systems, defined here as physical pairs 
with mass ratio in the range of 0.04 < M secondary I 'M primary < 0.2. 

This paper presents the results for the system AM 2306- 
721. This pair is composed by a peculiar spiral with disturbed 
arms (hereafter, AM 2306A) in interaction with an irregular galaxy 
(hereafter, AM2306B). Both galaxies contain very luminous Hn 
regions with Ha luminosity in the range of 8.30 x 10 39 <L(Ffa) < 
1.27 x 10 42 erg s -1 as estimated from Ho- images dFerreiro et al.l 
l2008h: and high star formation rate in the range of 0.07 to 10 M©/yr 
dFerreiro et alj2008h . The present paper is organized as follows: in 
Section(2] we summarize the observations and data reduction. Sec- 
tion [3] describes the gas kinematics of each galaxy and Section [4] 
present the numerical N-body/hydrodynamical simulations of the 
interaction. Section[5]presents the stellar population synthesis. The 
metallicity analysis is presented in Section [6] and the conclusions 
are summarized in Section|7] 



2 OBSERVATIONS AND DATA REDUCTION 

Long slit spectroscopic data was obtained in 20/21 June 2006 and 
23/24 June 2007 with the Gemini Multi-Object Spectrograph at 
Gemini South, as part of poor weather programmes GS-2006A- 
DD-6 and GS-2007A-Q-76. Spectra in the range 3 350 to 7 130A 
were taken in two settings with the B600 grating, and the 1" slit, 
yielding a spectral resolution of 5.5A. The frames were binned on- 
chip by 4 and 2 pixels in the spatial and spectra directions, respec- 




Figure 1. GMOS-S r'-band image of AM 2306-72 1 with the observed slit 
positions. 



Table 1. Journal of observations 



Date (UT) 


Exposure time(s) 


PA(°) 


AA(A) 


2006/06/20 


4x700 


190 


4280-7130 


2006/06/20 


4x700 


238 


4280-7130 


2007/06/24 


4x700 


190 


3350-6090 


2007/06/24 


4x700 


238 


3350-6090 


2007/06/24 


4x700 


118 


3350-6090 



tively, resulting in a spatial scale of 0.288 "pxl~\ and 0.9Apxl~' 
dispersion. 

Spectra were taken at three different position angles on the 
sky: PA=238° and PA=118° corresponding to positions along 
the major axis of AM2306A and AM2306B, respectively; and 
PA=190°, a position cutting across the disk of both galaxies. The 
exposure time on each single frame was limited to 700 seconds to 
minimize the effects of cosmic rays, with four frames being ob- 
tained for each slit position to achieve suitable signal. The slit posi- 
tions are shown in Fig.Q] superimposed on the r'-band image of the 
pair. Table [T] gives the journal of observations. Conditions during 
both runs were not photometric, with thin cirrus and image quality 
in the range 1.0" to 2.5" (as measured from stars in the acquisition 
images taken just prior to the spectroscopic observations). 

The spectroscopic data reduction was carried out using the 
gemini. gmos package as well as generic IRAF tasks. We followed 
the standard procedure for bias correction, flat-fielding, cosmic ray 
cleaning, sky subtraction, wavelength and relative flux calibrations. 
In order to increase the signal-to-noise ratio, the spectra were then 
extracted summing over six rows. Each spectrum thus represents an 
aperture of 1" x 1.73". For a distance of 119Mpc for AM2306A, 
and 116Mpc for AM2306B, (estimated from the radial velocities 
derived in Section[3]and adopting H =75 km s -1 Mpc~'), this aper- 
ture corresponds to a region of 577 x 998 pc 2 for AM 2306A; and 
562 x 972 pc 2 for AM 2306B. The nominal centre of each slit was 
chosen to be the continuum peak at A 5525 A. Fig. [3] and [4] shows 
a sample of spectra of AM 2306A and AM 2306B extracted along 
the slit for PA=238° and PA=1 18°, respectively. 




AM2306A 




Figure 2. Ha image of AM 2306B (left) and AM 2306A (right). Description and detailed analysis of these images have been presented in FP04. 



3 IONIZED GAS KINEMATICS 

The radial velocity was estimated from the strongest emission lines 
present in the spectra, namely Hy ,14340, Hp ,14861, [Om] ,15007, 
Ha ,16563, [Nn] ,16584, and [Sn] ,16717. The final radial veloc- 
ity for each spectrum was obtained by averaging the individual 
measurements from the detected emission lines, and the errors es- 
timated from the standard deviation of these measurements around 
that mean. 

Using the acquisition images in the r' filter, we did a simple 
isophote fitting using the stsdas. ellipse task to obtain the position 
angle of the major and minor axes of both galaxies. We found that 
the major axis for AM2306A is at PA=236°, and for AM2306B at 
PA=108°. The inclination of each galaxy with respect to the plane 
of the sky was also computed as cos(/) = bja, where a and b are the 
minor and major semi-axes of the galaxy, respectively. We obtained 
i = 56° for the main galaxy and i = 60° for the secondary, in 
agreement with the values measured by FP04. 

The rotation curves along the observed slit positions for both 
galaxies are presented in Figure(5] as well as the coresponding two- 
dimentsional spectra of Ha and/or H/3 emission lines. The radial 
velocity values after subtraction of the systemic velocities (as de- 
termined from the rotation model described below) are listed in 
Tabled 

Both galaxies show a fairly symmetric rotation curve, so we 
adopted a very simple approximation for the observed velocity dis- 
tribution, assuming that the gas moves under a logarithmic gravita- 
tional potential, following circular orbits close to a plane P(i, if/ ), 
characterized by its inclination to the plane of the sky (i) and the 
position angle (PA) of the line of nodes tp . This assumption results 
in an obs erved radial c i rcular velocity v(r, if/) in the plane of the sky 
given bv lBertola etal] dl99lh : 

, w V i?cos(i/r-i/'o)sin(Ocos(0 

V(r, f) = V s + ====== , (1) 

^R 2 rj + R 2 cos 2 (0 

with 



Table 2. Kinematical parameters 



Parameter 


AM2306A(PA=238°) 


AM2306B (PA=118°) 


Model 


Model 


V s (km/s) 


8919±5 


8669±3 


V (km/s) 


175±6 


185 ±6 


R c (kpc) 


[±0.3 


1.6 ±0.2 



r\ = [sin 2 (i/r - (fr ) + cos 2 (0 cos 2 (i/r - f )]> (2) 

where V, is the systemic velocity, R is the radius in the plane of the 
galaxy and V and R c are parameters that define the amplitude and 
shape of the curve. The fit of the rotation curves for each galaxy are 
shown in Fig.[5]and the parameters obtained are listed in Table|2] 

The heliocentric velocity of the main galaxy is found to be 
8919 km s -1 . The rotation curve obtained for AM2306A along 
PA=238° is well represented by the simple model above. This rota- 
tion curve is typical of spiral disks, rising shallowly and flattening 
at an observed amplitude of 145 kms~'. As can be seen in Fig(5] 
there are pronounced velocity features indicating local deviations 
from a smooth rotational field. One of such well defined features is 
located towards the northeast along PA=238°, at about 8 kpc pro- 
jected distance from the nucleus, with an amplitude (peak-to-peak) 
of about 100 kms -1 . This feature is located where the slit crosses 
the the outer spiral arm as seen in the Ha image (Fig. [2] right). This 
image also makes more obvious the disturbed spiral structure of the 
main component, very likely due to the interaction with AM 2306B. 

For the secondary galaxy, the rotation curve model results in 
an heliocentric velocity of 8 669 km s . The observed radial veloc- 
ities along PA=118° are also very well represented by the model. 
The rotation curve is quite similar to that of AM 2306A, rising shal- 
lowly and flattening at an observed amplitude of about 158 km 
s -1 . Although the young star forming population of the secondary 
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Figure 3. Spatial profiles of H/3, [Ora] ,15007 and ,15525 flux along PA=238° for AM2306A (left); and a sample spectra in the range of 3600 to 6000 A from 
different regions along PA=238° for AM2306A as marked in the shaded areas in the spatial profiles (right). The spectra have been normalized at ,15870 A 
and are plotted in the rest frame wavelength. 

Table 3. Radial velocities. The full table is available as Supplementary Material to the online version of this article from http://www.blackwell-synergy.com 



AM 2306A 






AM2306B 




PA=238" PA= 


190° 


PA= 


118° PA= 


190" 


R (kpc) V (km/s) R (kpc) 


V (km/s) 


R (kpc) 


V (km/s) R (kpc) 


V (km/s) 



13.96 SW 204±10 
12.96 SW 131±8 

11.96 SW 150±5 

10.97 SW 118±16 
9.97 SW 134±16 



13.96 S 
12.96 S 

11.96 S 

10.97 S 

9.97 S 125±U 



galaxy show clear morphological irregularities (Fig. [2] left), the ro- 
tation curve does not present significant deviations from the smooth 
repr esentation of the veloc i ty fiel d. 

iDonzelli & Pastorizal i 19971) estimated systemic velocities of 
8 762 and 9 069 km s" 1 for AM2306A and AM2306B, respectively. 
Following our results, we verified that the values presented by those 



authors were reversed, and therefore our estimates are found to 
agree to within 2 % with those previous determinations. 

An estimate of the dynamical mass can be derived by assum- 
ing that the mass inside a certain radius is given by M(R) = RV 2 /G. 
For the main galaxy, the deprojected velocity amplitude is 175 km 
s and its dynamical mass is 1.29 x 10 n Af© within a radius of 18 
kpc. For the companion galaxy, the deprojected velocity amplitude 
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Figure 4. Same as Fig. [3] but along PA=118° for AM2306B 



is 185 km s~' and the dynamical mass within a radius of 10.7 kpc 
is M(R) = 8.56 x 10 10 Mq. It is important to emphasize that the 
maximum radius to which we can observe the gas in emission is 
quite certainly smaller than the total radius of the galaxies, so our 
estimates of the dynamical mass give a lower limit to the actual 
dynamical mass of each system. 

In determining the dynamical mass, the main source of error 
is in the estimation of the deprojected velocity, which in turn is 
highly dependent on the assumed inclination of the galaxies with 
respect to the plane of the sky. For AM 2306A the inclination an- 
gle is well determined from the isophote fitting, with an estimated 
error of about 1%, propagating to a 2% uncertainty in the resulting 
total mass. For AM 2306B, the irregular morphology of the disk in 
optical images implies that the calculated value of the inclination 
angle varies significantly with radius. In this case, if we considered 
the isophote fitting from the regions more perturbed of AM 2306B 
to estimate the inclination angle, a maximum angle inclination of 
about 12 % would be obtained, implying in a reduction of up to 18 
% in the total mass. 



4 NUMERICAL SIMULATIONS 

In order to reconstruct the history of the AM 2306-721 system and 
to predict the evolution of the encounter, we modeled the interac- 
tion between AM2306A and AM2306B through numerical sim- 
ul ations using the the N-body/SPH code GADGET-2 developed 
bv lSpringel|_ (2005). T he gal axies were modeled following the pre- 
scription of Hernquistl ( ll993 l), including a gaseous disk component. 
The observed morphology and rotation curves presented in Sec. [3] 
provide the constraints to the simulations. 

As for any study of this kind, recreating the evolution of the 
AM 2306-721 system requires solving the reverse problem of cal- 
culating the orbit followed by the two galaxies from their observed 
properties. This, like most similar cases, is not a fully determined 
problem, since the observations do not provide all the necessary 
information to uniquely identify the solution. Therefore, in order 
to restrict the parameter space when setting up the initial condi- 
tions for the simulations, we first calculate orbits that satisfy the 
constraint given by the observed radial velocity difference, testing 
different eccentricities, pericenter and line-of-sight direction dis- 
tances. From that subset, based on the observed morphology, we 
select a few orbits to run the full simulations, from which the one 
that best fits the observed properties is selected. 

The orbit that best reproduces the observational properties is 
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Figure 5. Observed mean radial velocity and 2D spectra of Ha and Hfi emission lines as a function of apparent galactocentric distance. The left panels 
show the kinematics along PA=238°(top) and PA= 190° (bottom) of AM2306A, respectively. The right panels show the kinematics along PA=1 18°(top), and 
PA= 190° (bottom) of AM2306B, respectively. The velocity scale correspond to the observed and model values after subtraction of the systemic velocity of 
each galaxy, without correction by the inclination in the plane of the sky. 



found to be hyperbolic, with an eccentricity e = 1.15 and peri- 
galacticum of q = 5.25 kpc . The orbital plane is perpendicular to 
the plane of the sky, and intersects the later in the line that con- 
nects AM 2306A and AM 2306B as projected in the sky. In peri- 
galacticum, AM 2306B was behind AM 2306A, along the line of 
sight. The parameters for the best fit model are presented in Ta- 
ble|4] The composite rotation curve, as well as the individual model 
components and observed circular velocity, corrected by the incli- 
nation of each galaxy as given in Section[3] are shown in Figure[6] 
The models indicate a total mass of 1 .82 x 10 1 1 M© for AM 2306A 
and of 8.68 x 1O 1O M for AM 2306B. The mass of each individual 
component (disk, bulge, halo and gas) is quoted in Table|4] 

Figure [7] shows the time evolution of the encounter. Simula- 
tion starts 360 Myr before perigalacticum, and the time is shown 
in Myr in the upper right corner of each frame. The situation that 
best reproduces the morphology and kinematics of AM 2306-721 
system at the current stage is t = 610 Myr, or about 250 Myr after 
perigalacticum. The overall large scale morphology and kinematics 
agree well with observations, considering the low resolution of the 
simulation. 



5 STELLAR POPULATION SYNTHESIS 

A detailed study of the star formation in minor merger galaxies is an 
important source of information not only on the age distribution of 
their stellar population components, but to better understand several 
aspects related to the interacting process, its effect in the properties 
of the individual galaxies and their later evolution. The absorption 
features arising from the stellar component also affect to different 
degrees the measured intensity of the emission line in the spectrum 
of the gaseous component. This effect is more prominent in, but not 
restricted to, the Balmer lines, so the stellar population contribution 
must be subtracted from each spectra in order to investigate the 
physical properties of the gas in these galaxies. 

To investigate the star formation history of AM2306A and 
AM 230 6B we use t he stellar population synthesis method devel- 
oped by iBical dl998h . This method employs the equivalent widths 
W A of several spectral absorption features and the measured contin- 
uum fluxes F A at different wavelengths, comparing them to those of 
a model computed from a base of simple stellar population (SSP) 
elements with known ages and metallicities. The algorithm gener- 
ates all possible combinations of the base elements according to a 
given flux contribution step and compares the resulting W,\ and con- 
tinuum points to the input ones. The allowed solutions are those 
which reproduce, within predefined limits, the observed Wa and 
Fa- All these possible solutions are then averaged, and this aver- 
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Figure 6. Circular velocity curves for galaxy models AM 2306A (left) and AM 2306B (right). The rotation curves of the individual model components are 
shown: halo in orange, bulge in red, stellar disk in blue and gas disk in yellow. The composite rotation curve is the continuous thick black line. Points with 
error bars connected by thin black lines are the observed radial velocity curves, folded around the rotational centre and corrected by the inclination of each 
galaxy. 




Figure 7. Four snapshots representing the evolution of the encounter. Time is shown in Myr in the upper right corner of each frame, with respect to the 
beginning of the simulation. Stars are plotted in red, gas in blue. Orbital perigalacticum occurs at t = 360 Myr (upper right frame). Best fit to current 
evolutionary stage of the AM 2306-721 system is about 250 Myr after perigalacticum (bottom right frame). 



age is adopted as the final synthesis, with the uncertainty for each 
individual age component given by the standard deviation of the 
contribution from each allowed solution around the mean adopted 
value. The code used here is ba s ed on an upgraded version of the 
one presented bv lSchmittetai] Jl996l) . and includes a correction 
for the effect of the internal extinction. 

The spectral windows for measuring Wa are those defined i n 
the Lick system dWorthev & Ottavianilll997h iTrager et al.l 1998). 
with the addition of the K Ca n (corresponding to the spectral win- 
dow A A 3925 - 3945A)), and HI 1 A 3770A lines (spectral window 
A A 3765 - 3785 A). Previous to the measurements of the W A and F A 
the spectra were corrected by the f oreground (Galactic) reddening 
of E(B - V) = 0.03 mag taken from lSchlegel et al.l dl998h : and nor- 
malized to A 5870 A. The Wa were me asured using the PACCE code 
kindly provided bv lVale et alj J2007I) . The W A values measured for 
each galaxy at different galactocentric distances and position angles 
are listed in Table© 

The SSP base was created from the lBruzual & Charlotl f2003h 
evolutionary stellar population models, which are based on a high 
resolution library of observed stellar spectra. This library allows 
us to derive detailed spectral evolution of simple stellar popula- 
tions across the wavelength range of 3 200 to 9 500 A with a wide 
range of m etallicities. We used the P adova 1994 tracks as recom- 
mended by iBruzual & Charlotl J2003b . with the initial mass func- 
tion of Salpeter l Salpeteill 19551) . The final base contains five spec- 



tra corresponding to SSPs with ages of 2.5 Myr, 200 Myr, 1 Gyr, 
5 Gyr and 10 Gyr, and solar metallicity. The adopted number of 
elements in the base was defined by a compromise between hav- 
ing enough age resolution, having a limited number of features 
available to constrain the fit, and the desire to mantain a consis- 
tent base for all apertures in both galaxies. For the few positions 
with higher S/N spectra in each galaxy, we tested a more detailed 
age grid, including components with ages 2.5 Myr, 100 Myr, 200 
Myr, 500 Myr 1 Gyr, 5 Gyr and 10 Gyr, and solar metallicity. All 
those converged to the same results obtained with the original five 
component base, with the contributions from the 100 and 500 Myr 
populations resulting always less than 5 %. 

Figure [8] show an example of the observed spectra corrected 
for the reddening, the synthesized spectra and the pure emission 
spectra. Table ©presents the result of the stellar population synthe- 
sis for the individual spatial bins in each galaxy, stated as the per- 
centual contribution of each base element to the flux at /15870A. 
The values of E(B - V) derived from the synthesis are also given in 
Table© 

The spatial variations in the contributions of the stellar pop- 
ulation components are shown in Fig. [9] Large variations in the 
contribution from the different age components can be seen across 
the disk of both galaxies. The central regions of AM 2306A are 
dominated by the old (5-10 Gyr) population, with some significant 
contribution from a young 200 Myr and intermediate 1 Gyr compo- 



Table 5. Equivalent widths of the absorption lines. The full table is available as Supplementary Material to the online version of this article from 
http://www.blackwell-synergy.com 









Equivalent width (A) 


R(kpc) 


HI 1/1 3770 


KCall 


G/14300 Fe i A 5270 Fe 1/15709 Fe 1/15782 NaD/15890 


AM2306A(PA=238°) 



3.99 SW 
2.99 SW 
1.99 SW 
1.00 SW 




2.85±0.40 
2.32±0.19 



U2±0.12 



5.46±0.09 ... 2.35±0.03 

6.66±0.07 3.80±0.20 2.41±0.07 

5.69±0.23 4.00±0.31 2.15±0.08 

4.62±0.09 3.10±0.28 1.99±0.01 



0.69±0.05 1.47±0.01 

0.60±0.03 1.50±0.01 

0.61±0.02 2.15±0.02 

0.72±0.04 3.31±0.02 

0.75±0.02 4.29±0.02 




^ 0.5 - 



4500 5000 

Wavelength (A) 



C 2 



^ 1 



o ijlf 




i I 

4500 



5000 
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Figure 8. Stellar population synthesis for central bin along the major axis of AM 2306A (left) and AM 2306B (right). Top panel: spectrum corrected for 
reddening and the synthesized spectrum (shifted by a constant). Bottom panel: pure emission spectrum. 



nent of the spiral arms and the disk . On the other hand, the stellar 
population in AM 2306B is overall much younger, being dominated 
by the 2.5 Myr, 200 Myr and 1 Gyr components, which are quite 
widely spread over the whole disk. The star formation episode oc- 
curred about 200 Myr ago in both galaxies could be related with 
the perigalactic passage. 



6 EMISSION LINE INTENSITIES AND O/H 
ABUNDANCE 

Once the stellar population contribution has been determined, the 
underlying absorption line spectrum can be subtracted to allow 
the correct measurement and analysis of the line emission from 
the gaseous component. The line intensities were measured using 
Gaussian line profile fitting on the pure emission spectra. We used 
the IRAF splot routine to fit the lines, with the associated error 
being given as a 2 = cr 2 com + cr 2 ine , where cr cmt and cr nm are the 
continuum rms and the Poisson error of the line flux, respectively. 



The residual extinction associated with the gaseous component for 
each spatial bin was calculated comparing th e observed Hy/H/3 
and Ho-/H/? ratios to the theoretical values in iHummer & Storevl 
( 1987) for an electron temperature of 10 000 K and a density of 100 
cm 4 . The observed emission line intensities were then corrected by 
this residual extinction using the lHowarthUl983l) reddening func- 
tion. Table[6]present the main emission line intensities normalized 
to the flux in the H/J line. We derived the electron densities for 
AM2306A from the [S n]A 6717//1 6731 intensity ratios, using the 
temden routine of the nebular package from STSDAS/IRAF, assum- 
ing an electron temperature of 10000 K. The energy levels, tran- 
sition probabilities and collisional stre n gth values for [S n] were 
taken respectively from iBo wen ( 1960), K eenan et al.l dl993f) and 
iRamsbottom et alj dl996r) . The resulting electron density obtained 
for AM 2306A is in the low density limit for the [S n] lines, with 
values b elow N e = 630 cm 3 , compatible with the results for both 
,alactic iCopetti et al.l l200ol) . and giant extragalactic Hn regions 
CastanedaetalJl992h . 

The spatial metallicity variations across the disk of spiral 



Table 6. Stellar population synthesis results 



Flux fraction at A 5870 A 



R(kpc) 


2.5 Myr 


200 Myr 


1 Gyr 


5 Gyr 


10 Gyr 


E(B - V) 






AM2306-721A(PA= 


=238°) 






3.99 SW 


15±1 


52±2 


3+12 


15±62 


15±43 


0.42±0.21 


2.99 SW 


2±1 


14±1 


77±9 


1±I 


6±10 


0.16±0.03 


1.99 SW 


2±1 


1±1 


76±14 


13±42 


8±25 


0.13±0.02 


1.00 SW 


2±1 


8±1 


29±4 


16±46 


45±28 


0.18±0.05 





2±1 


22+1 


7±9 


13±72 


56±47 


0.29±0.10 


1.00 NE 


1±1 


23±3 


1±1 


5±27 


70±24 


0.06±0.01 


1.99 NE 


3+1 


1±2 


79±8 


12±22 


5+13 


0.14±0.03 


2.99 NE 


3±1 


21+2 


75+1 


1±1 


0±1 


0.20±0.05 


3.99 NE 


2±1 


24±6 


72±9 


1±2 


1±1 


0.07±0.01 


4.98 NE 


Oil 


30±1 


49+11 


7±19 


14±2 


0.02±0.01 


5.98 NE 


0±1 


39±2 


54±18 


4+11 


3±7 


0.13±0.03 


6.98 NE 


0±1 


35±6 


63±10 


1±2 


1±1 


0.14±0.03 


7.98 NE 


Oil 


50±1 


31 + 10 


9±32 


10±27 


0.22±0.06 






AM 2306-721B (PA= 


=118°) 






2.91 SE 


32±3 


67±6 


1±1 


0±0 


0±0 


0.40±0.46 


1.94 SE 


43±3 


55+5 


0±1 


1±1 


1±1 


0.35±0.37 


0.97 SE 


22±1 


53±2 


11+8 


7+18 


7+11 


0.39±0.18 





16+2 


47±4 


13±14 


14±54 


10±30 


0.23±0.07 


0.97 NW 


21±1 


45±2 


30±6 


2±5 


2±3 


0.39±0.18 


1.94 NW 


20±1 


41±2 


36±6 


2±5 


1+3 


0.31±0.12 


2.91 NW 


19±1 


42±2 


37±3 


1±3 


1+2 


0.28±0.09 


3.88 NW 


15±1 


44±2 


36±4 


3±4 


2±3 


0.16±0.04 


4.86 NW 


35±1 


63±3 


1±1 


1±1 


0±1 


0.31±0.12 


5.83 NW 


19±1 


44±1 


34±6 


2±5 


1±4 


0.01±0.01 


6.80 NW 


23+1 


25+1 


47±20 


3+11 


2±4 


0.02±0.01 


7.77 NW 


16±1 


46±2 


36±5 


1±4 


1±2 


0.01±0.01 



Table 7. Dereddened line fluxes 1(A). The full table is available as Supplementary Material to the online version of this article from 
http://www.blackwell-synergy.com 



AM2306A (PA=238°) 



R (kpc) 


c(Hy3) [On] ^3727 


H i 44340 


Hi i486 1 


[Ora] 44949 


[O m] 45007 


[Nn] 46548 


Hi 46562 


[Nn] 46584 


[Sn] 46717 


[Sn] 46731 


12.96 SW 







100±3 


40±2 


118±4 


22±4 


286±19 


60±5 


73±10 


73+10 


11.96 SW 


0.10 


47±9 


100±4 




72±4 


29±7 


286±25 


94±10 


94±14 


86+14 


10.97 SW 







100±5 




62±4 


54+13 


286±33 


117+17 






9.97 SW 







100±4 




55±3 




286±27 


121+15 






8.97 SW 







100±3 




61+3 


34±5 


286±18 


101+8 


89±4 


65±4 



Notes: all fluxes normalized to F(Hfi) as listed in Column 5. 



galaxies is a fundamental parameter for understanding galaxy 
evolution. In isolated spiral galaxies, it is often found that the 
metallicity in the ionized inter stellar medium decreases outwards 
dVila-Costas & Edmunds! 19921) . while a weaker or absent metallic- 
ity gradient in interacting galaxies of similar morphological type 
have been pointed out as evidence of radial gas inflows driven by 
tidal torques generated d uring the interaction dHenrv & Worthevl 
ll999l : lRov&Walshlll997h . 

Determining accurate element abundances from optical spec- 



tra is critically dependent on measuring temperature sensitive line 
ratios, such as [Om](/i4959 + A 5001) /A 4363. However, when 
doing spectroscopy of Hn regions with high metallicity and/or 
low excitation, temperature sensitive lines such as [Om]/i4363 
are found to be weak or unobservable, and empirical indicators 
based on more easily measured line ratios have to be used to es- 
timate metal abundances. The line ratio ^ 23=([On]/t3727+[ Om] 
A 4959+ [O m] 4 5007)/H/3, introduced by IPagel et all dl979h has 
been widely used for determination of oxygen abundance in 
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Figure 9. Synthesis results in flux fractions as a function of distance to the center along the major axis of AM2306A (upper panel) and AM2306B (bottom 
panel). 



Hn regions dVila-Costas & Edmund j 1 1992k IZaritskv et ail 1 19941 ; 
iGil de Paz et alj|2007h . 

Although the R23 indicator is sensitive to abundance vari- 
ations, it is double valued, presenting high to low metallic- 
ity turnover in the ran ge of 8 < 12 + log(0/H) < 8.3 
dPilvugin & Thuanll2005l) . Another critical issue regarding the use 
of the R23 is the fact that it is also sensitive to the ionization level 
of the emission gas, especially at low metallicities. These prob- 
lems can be solved, by choosing either the low or high-metallicity 
branch of R23, and combining it with another line indicator sen- 
sitive to the ionizati on level, such a s the ratio [On] (A 3727 + 
A 3729)/[0 m]^ 5007 jMcGaugbl[l99lh . 

Since we do not detect any temperature-sensitive emission 
lines in our spectra, we used the R23 indicator to estimate the metal- 
licity comparing the observed values with a gri d of photoioni zation 
models obtained from the code Cloudy/96.03 jFerlandll2002l). The 

[ rid w as built following the same procedures as iDors & Copettil 
2006h , with metallicities of Z = 2.0, 1.1, 1.0, 0.6, and 0.4 Z s , and 
ionization parameter log U = -2, -2.5, and -3. The solar va lue of 
12+log(0/H) = 8.69 is taken from lAllende Prieto et al.U200lh . For 
each model, the ionizing source was assumed to be a stellar cluster 
with energy distribution obtai ned using the s tellar population syn- 
thesis code 5 TARBURS T99 fceitherer et"al]|l999h . with an upper 



mass limit of 100 M B and an age of 2.5 Myr. Figure [10] shows the 
R23 versus [O n]/[0 111] diagramme, with the observed values super- 
posed in the computed models. Open squares correspond to regions 
in AM 2306A, filled triangles to regions in AM 2306B. Most of the 
points indicate gas with solar or less than solar abundance, with the 
exception of two spatial bins in AM 2306A, corresponding to the 
nuclear and the 1.0 kpc NE regions. Figures [TT1 and [T2l shows the 
estimated O/H abundance distribution as a function of the galacto- 
centric distance, R G , for AM 2306A and AM 2306B, respectively. 

The disk of the main galaxy (Fig. lilt shows a clear radial 
oxygen abundance gradient, that can be fitted as a linear func- 
tion 12 + log(0/H) = 8.75(±0.06) - 0.025(±0.007)R G , where 
Rq is in kpc. This abu ndance gradient is typical of spiral galax - 
ies, as found in M 101 dKinkel & Rosa! ! 1994- iGarnett et al.lll999h . 
and in NGC 1365 JPilvugin et al.ll2004l : lDors & C opetti 2005]). For 



AM 2306B the oxygen abundance is relatively homogeneous across 
the galaxy disk, not presenting a radial oxygen abundance gradient. 
The observed values are in the range 8.30 < 12 + log(0/H) < 8.58. 

The above results are in agreement with the luminosity- 
metallicity relation for gas rich galaxies (e.g. Zaritskv et alJI 19941 ; 



IContini et alj2002MLamareille et al.l2004l ; lTremonti et al.l2004h , in 
the sense that more luminous objects are more metallic. To investi- 
gate if the interaction affected the gas metallicity in both galaxies, 



Table 4. Parameters used on the simulations 





AM2306A 


AM2306B 


Number of points in disk 


2048 


2048 


Disk mass 


0.5 


0.17 


Disk radial scale length 


0.6 


0.6 


Disk vertical scale thickness 


0.06 


0.05 


Reference radius R re f 


0.06 


1.2 


Toomre Q at Rref 


1.5 


1.5 


Number of points in gas disk 


2048 


2048 


Gas disk mass 


U.Uj 


U.U1 / 


Gas disk radial scale length 


1.0 


0.75 


Gas disk vertical scale thickness 


0.02 


0.03 


Toomre Q at R,- e f 


1.5 


1.5 


Number of points in bulge 


256 


256 


Bulge mass 


0.08 


0.01 


Bulge radial scale length 


0.1 


0.06 


Number of points in spherical halo 


2048 


2048 


Halo mass 


3.5 


1.7 


Halo cutoff radius 


8.0 


4.0 


Halo core radius 


2.5 


0.6 



Notes: Simulations were done in a system of units with G=l. Model units 
scale to physical ones as: unit length is 3.5 kpc, unit velocity is 262 km s , 
unit mass is 5.586 X 10 10 M Q and unit time is 13.062 Myr. 

we compare the central O/H value of AM2306A a nd Am2 306B 
with i solated and pair galaxies from the sample of iKewlev et alj 
(2006). We found the central O/H v alue for AM2306A is in agree- 
ment with the results obtained by IKewlev et alj d2006l) for field 
galaxies at same luminosity, indicating that the interaction has not 
significantly affected the gas abundance in AM2306A. On the other 
hand, the central O/H value obtaine d for AM2306B is ab out 0.1- 
0.2 dex lower than those estimated bv lKewlev et alj 1 120061) for field 
galaxies at same luminosity, but is in good agreement with the val- 
ues obtained for galaxies in pairs. 

Interestingly, the mean value for AM 2306B, 12 + log(0/H) = 
8.39, is similar to the values found in the outer parts of AM 2306A. 
If we assume that in AM 2306B there was an oxygen abundance 
gradient before the encounter with AM 2306A, and that this gradi- 
ent was destroyed by gas flows from the outer parts to the center of 
the galaxy, mixing and h omogenizing the ch emical composition of 
the interstellar medium I Ke wlev et alj |2006). we can estimate the 
fraction of infalling gas that was required to produce the metallici- 
ties that are now observed in AM 2306B. For this, we assume that, 
before the interaction, AM2306B presented the same O/H gradi- 
ent observed now in AM 2306A. To reduce the initial abundance 
of 12 + log(0/H) ~ 8.75 in the central region to the current mean 
value of 8.39, it would be necessary that the central gas was diluted 
so as to contain a fraction of 56% of poor gas from the outer disk 
plus 44% of rich gas from the central region of the galaxy. 

Another possible explanation to mix the interstellar medium 
and flatten the r adial gradient would be the presence of a bar 
iFriedli et alj Jl994h . However, the very disturbed luminosity pro- 
file of this galaxy in B, V and / images (FP04) do not show any 
evidence of such feature. 

We can estimate the amount of ionized gas associated with the 
central star-forming com plex in AM 23 06B from the Ha luminosity 
of this region as given bv lTestorl fioOlh . For Rc, < 0.2kpc, we have 
L(Ha) = 2.71 x 10 41 erg s~' (FP04), and assuming a a filling factor 
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Figure 10. The relation log(S23) vs. log([0 n]/[0 m]) for the individual spa- 
tial bins in AM 2306A (squares) and AM 2306B (triangles). The curves rep- 
resent the photoionization models described in the text (dashed lines corre- 
spond to different values of the ionization parameter, solid lines to different 
gas metallicity). 
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Figure 11. The O/H abundance distribution as a function of galactocentric 
radius for the main galaxy along of the major axis. The solid line represents 
the fit of the data. 

of / = 0.30, we find a mass of ionized gas of M H u = 8.3 x 10 s Mq. 
According to the above proposed scenario, 56% of this total mass 
(4.6 x 10 5 Mq) would be infalling gas from the outer parts of the 
galaxy. Since we do not see any evidence of ionized gas flows in 
the radial velocities, we have to assume that the process of gas in- 
fall has already stopped. Considering the result of the numerical 
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Figure 12. The O/H abundance distribution as a function of galactocentric 
radius for the secondary galaxy along of the major axis. 

simulations, we take a timescale of ~ 100 Myr for the infall pro- 
cess (the gas motions were triggered by the perigalactic passage 
250Myr ago, but are no longer observable), resulting in an infall 
rate of 0.01Mo/yr. This infall rate is considerably smaller than the 
average gas inflow rate o f TM^/yr to the c entral 1-2 kpc predicted 
by recent merger models dlono e t al. 2004). 



7 CONCLUSIONS 

An observational study about the effects of the interactions in the 
kinematics, stellar population and abundances of the galaxy pair 
AM 2306-721 is performed. The data consist of long-slit spectra in 
the wavelength range of 3 350 to 7 130A obtained with the Gemini 
Multi-Object Spectrograph at Gemini South. The main findings are 
the following: 

(i) Rotation curves of the main and companion galaxies with 
an deprojected velocity amplitude of 175 km s -1 and 185 km s -1 
respectively were obtained. An estimate of the dynamical mass was 
derived for each galaxy, using the deprojected velocity amplitude. 
For the main galaxy, its dynamical mass is 1.29 x 10" M© within 
a radius of 18 kpc; and for the companion galaxy, the estimated 
dynamical mass is M(R) = 8.56 x 1O'°M0 within a radius of 10.7 
kpc. 

(ii) In the main galaxy, radial velocity deviations from the disk 
rotation of about 100 km s -1 were detected, which are probably due 
to the interaction with the companion galaxy. 

(iii) In order to reconstruct the history of the AM 2306-721 sys- 
tem and to predict the evolution of the encounter, we modeled the 
interaction between AM2306A and AM2306B through numeri- 
cal N-body/hydrodynamical simulations. The orbit that best repro- 
duces the observational properties is found to be hyperbolic, with 
an eccentricity e = 1.15 and perigalacticum of q = 5.25 kpc; the 
current stage of the system would be about 250 Myr after peri- 
galacticum. 



(iv) The spatial variations of the stellar population components 
of the galaxies were analysed by fitting combinations of stellar 
population models of different ages (2.5 Myr, 200 Myr, 1 Gyr, 
5 Gyr and 10 Gyr) and solar metallicity. The central regions of 
the main galaxy are dominated by the old (5-10 Gyr) population, 
with some significant contribution from a young 200 Myr and in- 
termediate 1 Gyr component along the disk of the galaxy. On the 
other hand, the stellar population in the companion galaxy is over- 
all much younger, being dominated by the 2.5 Myr, 200 Myr and 
1 Gyr components, which are quite widely spread over the whole 
disk. 

(v) Oxygen abundance spatial profiles were obtained using a 
grid of photoionization models and the i?23=([Oii]/i3727+[Oiii] 
A 4959+ [O m] A5007)/H/3 line ratio. The disk of the main galaxy 
shows a clear radial oxygen abundance gradient, that can be 
fitted as a linear function 12 + log(0/H) = 8.75(±0.06) - 
0.025(+0.007)Rg; in the companion galaxy the oxygen abundance 
is relatively homogeneous across the galaxy disk, with the observed 
mean value of 12 + log(0/H) = 8.39 

(vi) The absence of an abundance gradient in the secondary 
galaxy is interpreted as it having been destroyed by gas flows from 
the outer parts to the center of the galaxy. 
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SUPPLEMENTARY MATERIAL 

The following supplementary material is available for this article: 
Table 3. Radial velocities. 

Table 5. Equivalent widths of the absorption lines. 
Table 7. Dereddened line fluxes 1(A). 

This material is available as part of the online article from: 
http://www.blackwell-synergy.com/ (this link will take you to the 
article abstract). 
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